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Forthcoming photometric redshift surveys should provide an accurate probe of the acoustic peak 
in the two-point galaxy correlation function, in the form of angular clustering of galaxies within 
a given shell in redshift space. We investigate the form of the anticipated signal, quantifying the 
distortions that arise due to projection effects, and in particular explore the validity of applying the 
Limber approximation. A single-integral prescription is presented, which provides an alternative to 
Limber's equation, and produces a significantly improved prediction in the regime of interest. 

The position of the acoustic peak within the angular correlation function relates to the angular 
diameter distance to the far side of the redshift bin. Thicker redshift bins therefore shift comoving 
features towards smaller angular scales. As a result, the value of the photometric redshift error 
acquires a greater significance, particularly at lower redshifts. In order to recover the dark energy 
equation of state to a level of 1%, we find the total redshift dispersion must be determined to within 
Aa z < 10" 3 , which may prove challenging to achieve in practice. 



I. INTRODUCTION 

The propagation distance of sound waves prior to de- 
coupling is a convenient reference scale embedded in the 
distribution of galaxies. Lately these Baryon Acous- 
tic Oscillations (BAO) have received particular atten- 
tion due to their potential to measure the dark energy 
equation of state, w 0, @, H, 0]. To attain competitive 
constraints on w, the angular scale of the BAO must be 
determined at the sub-percent level. In addition, robust 
models are required for any phenomena capable of in- 
ducing systematics of a comparable magnitude. In the 
near future, large photometric redshift surveys such as 
Pan-STARRS T| and DES will cover a substan- 
tial proportion of the sky, providing a vast map of large 
scale structure in the universe. However, the precision 
and reliability of the photometric redshifts will signif- 
icantly influence the dark energy constraints. With a 
typical rms redshift error of a z ~ 0.05 the longitudinal 
modes of the power spectrum are exponentially damped, 
leaving only those which are closely aligned to the plane 
of the sky. Therefore BAOs from photometric redshift 
surveys are predominantly sensitive to the angular di- 
ameter distance Da- This may also be interpreted as 
an effective reduction in volume given by approximately 
12(ct z /(1 + z)/0.03) [|, although Cai et al. @ reach a 
more optimistic conclusion, suggesting only a fivefold re- 
duction in volume. 

The characteristic signature for a particular dark en- 
ergy model resides in the apparent redshift variation of 
the angular acoustic scale. In order to study this be- 
haviour, the galaxy population may be divided into thin 
redshift slices, limited in thickness to be > a z . We ex- 
plore the smoothing effect which arises when projecting 
from three dimensions down to two. Locating the acous- 
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tic signature in Fourier space and real space - via the 
power spectrum and correlation function - are each sus- 
ceptible to different forms of systematic errors. Finding 
both in agreement would therefore strengthen our confi- 
dence in any subsequent constraints on cosmological pa- 
rameters. For this study we focus on the acoustic peak 
in the angular correlation function, paying particular at- 
tention to the influence of projection effects, and analyse 
the shift in peak position as a function of bin width. The 
Limber approximation, recently examined by Simon 0], 
acts as a starting point from which we develop an alterna- 
tive analytic description of the transformation from the 
spatial to angular correlation function. Recently LoVerde 
& Afshordi [8| express the angular power spectrum as 
a power series, within which the Limber approximation 
appears as the th order term. Our approach is some- 
what different, and we obtain a single- integral approxi- 
mation for the angular correlation function, which pro- 
vides a good description for galaxy redshift distributions 
of practical interest. 



II. FIDUCIAL COSMOLOGY 

To evaluate the spatial correlation function we be- 
gin with the linear matter power spectrum output from 
CAMB [9(, and adopt the default cosmological parame- 
ters Q m = 0.27; fl A = 0.73; fl b = 0.046; n s = 0.96; ct 8 = 
0.79. Throughout this work we assume a flat universe. 

Conventional means of modelling non-linear structure 
formation, such as the HALOFIT prescription [l(| are 
not well suited to an oscillatory P(k) since they focus on 
corrections in the highly non-linear regime, rather than 
the quasilinear region where the BAO signal resides. Re- 
cently Eisenstein et al. [|[ offered a simple approach to 
modelling the degradation of the acoustic peak. By con- 
sidering the rms radial displacement for a pair of galaxies 
separated by 100 h~ 1 Mpc, the damping is given by 






FIG. 1: A representation of the observable pairings in the case of (i) a thin redshift slice (ii) a thick redshift slice with truncation, 
as described in ([IV B[) and (iii) a thick redshift slice without truncation, as provided by the Limber approximation ([8]). Each 
point on the shell's surface - along with its corresponding antipode - represents a potential galaxy pair. 
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where D oc S(z)/S(0) is the growth function (normalised 
such that D = a at early times), and sq is taken to be 
10.9/i -1 Mpc, rescaled from the 12.4 ft, _1 Mpc quoted in 
Q to match our lower value of as . The real-space corre- 
lation function at a given redshift is then obtained by 
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However in redshift space the power is anisotropic, so 
before proceeding further we first explore the extent to 
which redshift distortions may impact our analysis. The 
anisotropic correlation function £(r_i_,r||), where r± and 
rii are the transverse and radial parts of r, is related to 
the power spectrum via 
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where [ik is the radial part of k. Thus e.g. £(rj_, 0) is not 
in general the same as the real-space correlation function. 
Incorporating a simple model of the Kaiser and Fingers 
of God effects leaves us with 
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where a v denotes the one-dimensional rms single-point 
velocity dispersion. The strong damping of radial modes 
ensures that for all but the smallest values of k, contribu- 
tions from the Kaiser effect are restricted to small values 
of [i. This resembles a slow-varying scale-dependent bias 
in the angle-averaged power spectrum, which is therefore 
unproblematic provided we adopt an empirical "wiggles 
only" fitting function. 

In practice we will want to consider angular cluster- 
ing in shells of width comparable to the photo-z error, 
o z ~ 0.05. This scale is large by comparison with those 
of peculiar velocities, so we shall ignore the above compli- 
cation and treat the situation as though photo-z's are a 
noisy measure of the true radius, even though this is not 
strictly true. We shall thus tend to use the term 'redshift 
shell' to mean 'radial shell'. 



III. PROJECTION EFFECTS 

To begin, we consider an infinitesimally thin radial 
shell, so clustering may only be viewed perpendicular to 
the line of sight. In this case the angular correlation 
function simply resembles a rescaled version of the spa- 
tial correlation function, Wt{0) ~ £(Xs^)j where Xs is the 
(comoving) radial distance to the bin. This is known as 
the thin-layer approximation. Any given part of the spa- 
tial correlation function could therefore be represented as 
a circle on the sky, as illustrated in the left hand panel 
of Fig. [U 

When considering a thicker redshift bin, with a uni- 
form selection function, two distinct phenomena influ- 
ence our observation. The first, as represented by a trun- 
cated sphere in the centre panel of Fig. [TJ is that galaxy 
pairs with significant line-of-sight displacements begin to 
contribute to the signal. Therefore any given large-scale 
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FIG. 2: The angular correlation function (solid) for a top-hat galaxy redshift distribution evaluated at z = 0.5, compared to the 
Limber approximation (dashed), and thin-layer approximation (dotted). The bin widths are scaled such that Az — 0.1(1 + z) 
(left) and Az = 0.04(1 + z) (right). 
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FIG. 3: The same plots as in Fig[2] but with the bin located at z = 1. 
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FIG. 4: The same plots as in Fig[2] but with the bin located at z = 3. 
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feature from £ (r) begins to influence the smaller angular 
scales of w(6). The extent of this is governed by the ratio 
of the clustering scale to the bin thickness. The second 
effect arises since the truncated sphere subtends a differ- 
ent angle depending on its location within the bin. This 
effect gains prominence when the bin thickness is a sig- 
nificant proportion of the radial distance to the bin. This 
can often be alleviated by adopting a projected correla- 
tion function rather than an angular correlation function, 
but the presence of photometric redshift error proves to 
be a limiting factor. 

For the case of a redshift bin much thicker than the 
clustering scale, as seen in the right hand panel of Fig. [TJ 
our model is simplified since any truncation can be ne- 
glected and only the projection effect persists. Therefore 
any galaxy pair whose centre lies within the bin is as- 
sumed to contribute to the signal. This simplification 
leads to the erroneous inclusion of pairs which extend 
beyond the edge of the bin. Indeed this is the situation 
described by the Limber approximation, which provides 
us with a good starting point for the purposes of study- 
ing the acoustic peak. The clustering scale of interest 
is around 150 Mpc, whilst photometric redshift precision 
ensures the effective bin thickness is typically double this. 



A. Limber Shift 



This simplification utilises the flat-sky approximation, 
which holds well even on the relatively large angular 
scales under consideration. At a separation of 5° the ap- 
proximation (1 + cos#) ~ 2 provides an overestimation 
of 0.2%, and this propagates to R with an error of less 
than 0.1%. This is more prominent for the longitudinal 
modes, to which we are least sensitive. 

Figs. [2HH illustrate the shift in the acoustic peak as 
predicted by the Limber equation. Note how the acous- 
tic peak is almost entirely erased when the bin width is 
Az — 0.1(1 + z), highlighting the need for accurate pho- 
tometric redshifts. The increasing bin thickness, chosen 
to mimic the degradation of the photometric redshifts, 
is counteracted by the reduced value of dx/dz at higher 
redshift, leading to little change in the projection effect 
across the range of redshifts 0.5 < z < 3. Ordinarily the 
Limber approximation overestimates the signal, however 
an underestimate arises here due to the influence of the 
correlation function's negative regime. The Limber ap- 
proximation clearly offers a more appropriate estimation 
than that of the thin layer. Its precision varies according 
to the nature of the chosen fitting algorithm, and im- 
proves with broader bins. We find the relevant angular 
scale is underestimated, with the magnitude of the dis- 
crepancy typically ranging from 2% to 5%, substantially 
exceeding the level of accuracy required for dark energy 
studies. 



The angular correlation function is evaluated numeri- 
cally using (23) from Simon Q 

l + cosfl JJ dy v -§dy y +§ x 

(6) 

where the variables R, y, and x denote the galaxy pair 
separation, their mean distance from the observer, and 
their line-of-sight separation, as given by 

/ 2R?~ V(l- cos gj 

x -y i+cos0 ' (n 

The evolution of £ is represented by its first argument, 
whilst the second provides the galaxy pair separation R. 
This formulation assumes that only linear evolution of £ 
occurs over the distance x, and no significant change in 
position of the nearer galaxy arises during the light-travel 
time between the pair. 

Projection effects lead to the acoustic peak appearing 
at a smaller angular scale than one might naively expect. 
Here we assess whether a simplified form of © may be 
used to provide an estimate of the acoustic peak location. 
The Limber approximation is given by 

/>oo poo /j \ 2 

w(9)= J J (—) t(y,Vx 2 +y 2 2 ) dxdy. (8) 



IV. SHIFT ESTIMATION 

As we have seen, adopting the Limber approximation 
results in a good but insufficiently accurate estimation for 
the positioning of the acoustic peak. Here we attempt to 
gain a quantitative understanding of the shift, and inter- 
pret this as a smearing of Wt(0) 1 the correlation function 
from an infinitesimal bin. 



A. Limber Equation 

We begin by transforming the Limber equation to 
plane polar coordinates 

W{9)= 8 L I {dy^J ^ rdrd6 '> ( 9 ) 

where we substitute x = rsm9', y9 = rcos8'. 

For the moment we neglect the redshift evolution in the 
correlation function, and assume a top-hat galaxy dis- 
tribution, simplifying our equation by assuring dn/dy = 
(A%) -1 within the top-hat, where A% corresponds to the 
physical depth of the redshift bin. Later we shall see that 
our results are largely insensitive to the functional form 
of the redshift bin, and the impact of an evolving £(x) 
across reasonably narrow redshift bins is minimal. 

These simplifications leave us with an angular correla- 
tion function given by 
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and it is instructive to study the case where the correla- 
tion function consists of a delta function at an arbitrary 
location ro : 

£(r) = S(r - r), (11) 



«*> (12) 

where the lower integration limit is given by 

0[ = arccos | ) , (13) 
V ro ) 



Xi=Xs + ^, (14) 

and similarly for 9' 2 : 
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FIG. 5: The solid line represents the smear function derived 
from the Limber approximation, w 3 as defined in (|18p . and 
evaluated at z — 0.5 with Az = 0.1. The dashed line is the 
exact expression, illustrating the modification introduced by 
accounting for the finite bin width. On smaller angular scales, 
corresponding to a large line-of-sight separation, the pairs lost 
from truncation become particularly apparent. Of the three 
vertical dotted lines, the central one represents the original 
<5-function, while the left and right lines denote the angular 
scales associated with the far and near bin edges respectively. 
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where we denote the upper and lower boundaries of the 
galaxy distribution by xi and Xi respectively, and define 
arccos(x) = for x > 1. This leaves us with 
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Substituting ro = Xs@o leaves us with the Limber smear 
function 



At greater angular scales contributions are rapidly lost, 
starting from those at the far edge of the bin, with the 
signal reaching zero at the angular scale corresponding 
to the near edge of the bin. Conversely, smaller angular 
scales gradually receive fewer contributions since they re- 
quire the pair orientation to be increasingly aligned with 
the line of sight. These are the pairs which prove partic- 
ularly susceptible to overestimation by the Limber equa- 
tion, since they readily protrude beyond the bin bound- 
ary. 

Now returning to (|18p . and extending our analysis to 
a general correlation function, the parameter 9q becomes 
a variable over which we integrate 
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The functional form of this is illustrated in Fig. and 
this highlights the mechanism by which the peak position 
is shifted to smaller angular scales, typically by around 
10%. The projected signal peaks at 9 p , the angle cor- 
responding to the far edge of the redshift bin, as given 
by 
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B. Correcting for Bin Width 

Here we modify the limits of our integration in order 
to account for the finite bin width. This aims to improve 
on the Limber approximation's estimate for the location 
of the acoustic peak. 

Towards the bin edges the spherical shell assumption 
breaks down - a galaxy cannot pair up with another 
galaxy which lies beyond the bin boundary. This trunca- 
tion is particularly problematic on scales which are large 
relative to the bin width (see Fig. [1]). By considering 
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the point at which one of the galaxy pair extends be- 
yond the bin boundary we shrink the integration range 
accordingly with our new limits given by 
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FIG. 6: Here we take the angular correlation function corre- 
sponding to a Gaussian redshift distribution at z — 0.5 with 
<j z — 0.02(1 + z), and superpose that of a uniform galaxy 
distribution (dot-dash), scaled to provide a matching rms. 
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sin 02 = — ^-cos - . (26) 
The angular correlation function is then recovered by 
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This expression is exact within the context of a uniform 
galaxy distribution and flat sky approximation. 

The influence of this adjustment is illustrated by the 
dashed line in Figf5j The Limber equation generally over- 
estimates the amplitude of the correlation function due 
to over-counting the number of pairs, and this truncated 
integral provides the necessary correction. Upon evalu- 
ating the integral, we recover the solid and dashed lines 
originally shown in Figf2] 

V. THE GALAXY DISTRIBUTION 

Up until now we have considered a uniform galaxy dis- 
tribution. Provided Ax/x is small, the functional form 
of n(z) has a limited impact on the appearance of the 
correlation function. We find that the w(0) produced by 
a Gaussian selection function may be closely mimicked 



by simply selecting a top-hat selection function, of width 
Az chosen with a matching rms, given by 

Az = Vl2a(n(z)). (29) 

Fig. [6] compares three estimates of the angular corre- 
lation function for a Gaussian distribution. We find that 
this top-hat approximation compares favourably to the 
conventional Limber and thin layer approaches, produc- 
ing a more faithful reproduction of the bump's shape and 
magnitude. Yet for higher precision and larger values of 
Az/z, inevitably one must return to the exact form given 
by ©. 

VI. CALIBRATION REQUIREMENTS 

In this section we clarify what constraints on the pho- 
tometric redshift errors are needed in order to avoid sig- 
nificant degradation or biasing of the dark energy con- 
straints from future surveys. 

A. Photometric Redshift Error 

When projecting from the 3D correlation function to 
the two-dimensional angular correlation function, not 
only does the smoothing process dampen the appearance 
of the peak, but its angular position shifts in relation 
to the thickness of the redshift bin. This modified peak 
location, previously illustrated in Figure [5l is associated 
with the angular diameter distance at the far edge of the 
redshift bin. It is via this mechanism that the magni- 
tude of the photometric redshift error controls the peak 
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B. Dark Energy Constraints 

The process of fitting the acoustic peak to the data is 
not straightforward. A strict fitting regime correspond- 
ing to the anticipated shape of £(r) from simulations 
would leave us vulnerable to systematic errors. Mean- 
while the adoption of a fitting function with many degrees 
of freedom would substantially degrade the precision of 
our measurement. To avoid introducing any dependence 
on the fitting method into our results, we simply consider 
the angular shift induced on a correlation function of the 
form £(r) = S(ro — r), where tq = 110 /i _1 Mpc. 

In Figure [7] we evaluate the impact of an uncertainty 
in a z at redshifts z = 0.3,0.5,1, with fixed bin width 
Az p = 0.1 and a z — 0.05(l + z). We incorporate an error 
Aov defined such that 



FIG. 7: Overestimating the photometric redshift error may 
lead to a substantial overestimate of the dark energy equa- 
tion of state, particularly at lower redshifts. The vertical line 
denotes the achievable limit in the near-term future, as given 
by Newman [ill ]. 



position. Any uncertainty in a z therefore influences our 
estimate of the dark energy equation of state. 

To quantify this effect we select a top-hat of width 
Az p = 0.1 in photometric redshift space, which is then 
convolved with a Gaussian of standard deviation a z to 
simulate the true redshift distribution. As demonstrated 
in the previous section, the resulting angular correlation 
function is readily approximated by simply adopting a 
pure top-hat selection function of width 



Az = jAzl + 12a 2 z , 



(30) 



and the redshift of interest is z\, the far edge of the bin 



z\ = z + 



Az 



(31) 



It is therefore crucial to acquire accurate and unbiased 
estimates of both the mean redshift z and the photomet- 
ric error a z . Note that from (f3"Tj) we find 



dzi 
da z 



6crz 
Az 



(32) 



thus if we take Az p = 2<r z , this leaves our measurement 
50% more sensitive to the value of a z compared to the 
mean redshift z. A selection function much wider than 
2er 2 is undesirable as it leads to excessive damping of the 
peak, and leaves us susceptible to the redshift evolution 
of a z . 

To proceed, we assume the mean redshift is perfectly 
known, and explore the impact of a z on our estimation 
of the dark energy equation of state. 



Az 



12(a z + A(i z 



(33) 



Our estimate of the equation of state, w, is then per- 
turbed by an amount relating to dDA/dzi, thus leading 
to the steeper slopes at lower redshift. Provided the esti- 
mator is unbiased, averaging over several bins could im- 
prove the situation. Conversely, bypassing the CMB by 
taking the ratios of angular scales in different bins would 
further aggravate the problem. 



C. Measuring a z 

The conventional method for estimating the photomet- 
ric redshift error is to use a calibration set of spectro- 
scopic redshifts covering a suitably representative sample 
of galaxies. Newman [111 ] provided an estimate of Aa z 
attainable from spectroscopic calibration samples, com- 
bined with information from their angular clustering, and 
found Aa z ~ 2 x 10~ 3 . This includes a modelling of the 
outlier population as a second Gaussian with a signifi- 
cantly larger standard deviation. 

Zhan & Knox [3 consider using the exponential 
damping of the longitudinal modes of the galaxy power 
spectrum to determine the ratio cr z /H, where H is the 
Hubble parameter at the redshift of interest. Applying a 
prior on H will then provide an estimate of a z . For ex- 
ample, for an LSST-like survey, they find a 1% prior on 
H leads to an uncertainty on a z of 1%. Ma et al. [l3|] ex- 
plore the impact of photometric redshift uncertainties on 
weak lensing surveys, which appear to have photometric 
calibration requirements at a similar level to that of the 
BAO. 

The non-Gaussianity of photometric redshift errors, 
along with their strong variation with redshift (see Fig 8 
of [a]), ensures that we require more degrees of freedom 
than simply a z , which further heightens the difficulty of 
reaching the required level of precision. Ma & Bernstein 
[T3 ] estimate this may be compensated for by enlarging 
the spectroscopic calibration sample by a factor of five. 
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VII. THE PROJECTED CORRELATION 
FUNCTION 

Here we justify the choice of an angular correlation 
function, used throughout this work, as opposed to the 
projected correlation function. This is defined such that 

w P (r±) = 2 [°° £(r) J—^. dr. (34) 

Jr ± V r ~ r ± 

To evaluate this quantity we require the precise loca- 
tion of the galaxy in the bin, but this is unknown since 
a z ~ Az. As a result, adopting this approach aggra- 
vates the smoothing of the correlation function. To illus- 
trate this we consider the range of transverse separations 
which are observed by a particular value of w(9) and 
WpifiXs)- F° r the angular correlation function this is sim- 
ply given by the limits xi@ and Xi®- However with the 
projected correlation function this window is extended 
to (0xsXi/X2) and {OxsXz/xi) respectively. This occurs 
when pairs thought to be located at xi are actually at 
X2, and vice versa. 

Conversely when studying a bin with a z < Az the 
projected correlation function is the more natural op- 
tion, and will be subject to projection effects in a manner 
closely analogous to that considered here. 

VIII. DISCUSSION 

We present a single-integral solution to the angular 
correlation function, by considering the special case of 



a uniform galaxy distribution. This provides a useful 
alternative to the traditional Limber and thin-layer ap- 
proximations, and offers a substantial improvement in 
accuracy when evaluating the appearance and location 
of the acoustic peak. 

This top-hat approximation can be utilised to derive an 
analytic estimate for the impact of photometric redshift 
errors on BAO surveys. In agreement with Zhan & Knox 
[l2rj , we find the value of o z ought to be known to better 
than ~ 10~ 3 if the acoustic peak is to provide significant 
and robust constraints on the dark energy equation of 
state, and this is particularly apparent for lower redshifts, 
z < 1. It also seems likely that the anticipated departure 
from Gaussianity of this error distribution must be well 
understood. 

We have not addressed the complexities associated 
with galaxy bias, since the same projection process ap- 
plies regardless of the original form of the spatial corre- 
lation function. However, within the context of a photo- 
metric redshift survey a well-defined model of £(r) will 
be required in order to reverse the smoothing process. 

Whilst future projects such as LSST and SNAP will 
look to higher redshifts, it is the near-term projects such 
as Pan-STARRS and DES which, despite their higher 
tolerance for statistical errors, may be more susceptible 
to the difficulties in mapping the low-redshift behaviour 
of <r z (z). 
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